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ABSTRACT 


fit 


A  continuous  semigroup  of  linear  contractions  (S(t))  t  >  0  on  a  Banach' 
space  X  can  be  generated  through  the  product  formula 


(1) 


S(t)x  -  lim  u(— )nx  Vx  e  X  Vt  >  0 


n+°° 


where  (u(h)  )j1>0  is  any  family  of  contractions  on  X  whose  right-derivative 
at  h  =  0  coincides  with  the  infinitesimal  generator  A  of  S(t),  that  is 


(2) 


Vx  e  D( A)  lim 
h+0 


x  -  U(h)x 
h 


=  Ax  (=  lim 
t+0 


S(t)x 


)  • 


Actually  formula  (1)  extends  to  the  case  when  S(t)  is  a  semigroup  of 
nonlinear  contractions  "generated"  by  an  m-accretive  operator  A  and 
condition  (2)  may  also  be  weakened  to 

I  -  U(h)  h+Q  . 

(3)  - - -  — -*  A  in  the  sense  of  graphs. 

We  look  here  at  the  same  formula  ( 1 )  when  the  regular  step-size  t/n  is 
replaced  by  a  variable  step-size,  namely 

(4)  lim  u(hnn)  u(h°n  )  ...  u(h”)x  =  S(t)x  Vx  e  D( A) 

n+“  N  N  -1 

where  )  hn  =  t  and  lim  (  max  hn)  =  0. 
u  n  i  v  n  1 ' 

1<i<N  n+«  14i<N 

Surprisingly,  it  turns  out  that  (4)  fails  under  assumption  (3);  we 
exhibit  a  counterexample.  However,  we  prove  that  (4)  holds  true  under  the 
stonger  assumption  (2). 

AMS  (MOS)  Subject  Classification:  47H20,  65J15,  35R15. 

Key  Words:  Product  formulas,  nonlinear  semigroups,  accretive  operators, 
approximation  scheme,  stability  and  convergence. 
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SIGNIFICANCE  AND  EXPLANATION 


Many  numerical  schemes  are  based  on  the  following  abstract  product 
formula:  given  an  evolution  equation  which  can  be  written  in  the  abstract 

form 


4^(t)  +  Au(  t ) 
dt 


0  on  (  0 ,  T ) 


u( 0 )  =  uq  given 


with  A  being  a  suitable  operator  on  a  functional  space,  then  under  adequate 
assumptions  the  solution  u  is  given  by  the  formula 


u(t) 


lim  u(— )nu 


n-*» 


0 


where  (u(h) )  h  >  0  is  a  family  of  mappings  whose  right-derivative  at  h  =  0 
coincides  with  A,  that  is 

Vx  e  0(A)  lim  x.  gm*  =  AX. 
h+0  h 

This  a  natural  extension  of  the  exponential  formula 

u ( t )  =  lim  (i  +  -  A)“n  uQ  . 
n+“ 


Above  product  formula  provides  an  extensively  used  tool  for  the  Lie-Trotter 
formulas  insuring  the  convergence  of  various  alternating  direction  methods  for 
evolution  equations. 

In  view  of  numerical  applications,  it  is  important  to  relax  the  restric¬ 
tion  that  the  time  step-size  be  regular  in  the  above  schemes  when,  for  instance 
coupled  with  path  control  technics  or  even  for  stability  requirements. 

The  goal  of  this  paper  is  to  decide  whether  the  Chernoff  formula  is  valid 
when  the  step-size  is  variable.  It  turns  out  that  under  the  usual  assumptions 
instability  and  lack  of  convergence  can  appear  as  soon  as  the  step  size  is 
slightly  perturbed.  We  give  here  an  example  in  that  direction.  However,  we 
prove  that  slightly. stronger  assumptions  can  prevent  these  pathological 
situations.  -  Vao' >■ 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


ABOUT  PRODUCT  FORMULAS  WITH  VARIABLE  STEP-SIZE 
Michel  Pierre*  and  Mounir  Rihani* 

1.  Introduction. 

Given  (s(t))t>Q  a  continuous  semi-group  of  linear  contractions  in  a 
Banach  space  X,  it  can  be  generated  via  the  following  exponential  formula: 

(1)  Vt  >  0;  Vx  e  X,  S(t)x  =  lim  (i  +  -  A)"nx 

n-M»  " 

where  -A  denotes  its  infinitesimal  generator  and  I  is  the  identity.  More 
generally,  one  can  replace  in  (1)  the  resolvents  (I  +  h  A)-1  by  a  family  of 
contractions  (u(h))^>Q  from  X  into  X  satisfying 

(2)  Vx  e  D( A)  lim  x  ~  ^(h)x  =  Ax  . 

h+0  h 

Thus,  it  has  been  proved  (see  [11],  [3])  that  (2)  implies 

(3)  Vt  >  0,  Vx  e  X,  lim  b(-)nx  =  S(t)x  . 

n+“  n 

Actually,  as  proved  in  [7],  [4],  the  condition  (2)  can  be  weakened  to  assuming 
that 

(4)  Vx  e  X,  VA  >  0,  lim  (l  +  ^I~T(h|)  )~1x  =  (1  +  \A)_1x 

h-K)  h 

which  is  equivalent  to  saying  that  the  graph  in  X  x  X  of  the  operator 
h  ^(i-U(h))  converges  to  the  graph  of  A  as  h  tend  to  0. 

The  implication  (4)  =>  (3)  turns  out  to  have  a  nonlinear  version  as 
proved  by  Brezis  -  Pazy  [2] .  Namely,  let  A  be  an  m-accretive  operator  in 
X,  that  is  A  C  X  x  X  and 


*  Department  of  Mathematics,  University  of  Nancy  I,  B.P.  239,  54502- 
Vandoeuvre,  Les-Nancy,  France. 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-004 1 


(5)  VA  >0,  (I  +  Aa)  is  a  contraction  everywhere  defined  on  X. 

Let  (s(t))t>Q  be  the  nonlinear  semigroup  generated  by  A  in  the  sense  of 
Crandall  -  Liggett  [5],  that  is 

(6)  Vt  >  0,  Vx  e  D(A),  S(  t )  x  =  lim  (l  +  —  a)~Dx  . 

-  v  n  J 

Then,  if  (u(h))hJ>Q  is  a  family  of  contractions  from  X  into  X  satisfying 
(4),  we  have  the  product  formula 

(7)  Vt  >  0,  Vx  e  D( A) ,  lim  uf— )nx  =  S(t)x  . 

n 

n-n» 

The  purpose  of  this  paper  is  to  look  at  the  formula  (7)  when  the  regular 

step  size  —  is  replace  by  a  variable  step-size.  In  particular,  under  the 
n 

same  assumptions  as  above  does  one  have 

(8)  lim  u(h^  )  u(h"  _t)  ...  u(h^)x  =  S(T)x 

n>°°  n  n 

when 

Nn 

\  hn  =  T,  lim  (  max  hn)  =  0? 

i=1  1  n>°°  1<i<N  1 

n 

We  give  here  a  negative  and  a  positive  answer  to  that  question: 

а)  given  A  an  m-accretive  operator  the  condition  (4)  does  not  imply 
(8)  as  proved  by  a  counterexample  described  in  paragraph  3. 

б)  if  A  is  a  single-valued  m-accretive  operator  and  if  U(h) 
satisfies  the  stronger  assumption  (2),  then  (8)  holds. 

Acknowledgements .  The  authors  thank  Mike  Crandall  for  several  helpful 


liscussions  on  that  subject. 


2,  The  positive  result « 

Let  X  be  a  Banach  space  with  its  norm  denoted  by  | • |  and  let  A  be  a 
single-valued,  accretive  operator  on  X,  that  is 


VX  >  0,  Vx,  x  e  D(A),  |x  -  x|  <  |x-x+X(Ax-Ax)|  . 


For  h  e  (0 ,d)  let  U(h):  D(  A )  — ►  D( A)  satisfy 


Vh  e  (0,hQ),  Vx, x  e  D( A)  | U( h) x  -  U(h)x |  <  | x  —  x | 


Vx  e  D(  A) ,  lim  ^  —  =  Ax  . 

h+0  h 


Given  T  >  0,  0  a  subdivision  of  (0,T): 


(12)  0  tg  <  t^  <  •••  <  1  ^  T  ^  t«,  |o|  “  max  (t,  —  t.  .  ) 

1<i<N 

and  z:  [0,T[  +  X  such  that 


vi  «  i,...,n,  vt  e  jtj..,,  t±] ,  z(t)  *  zL  e  x, 


let  us  assume  there  exist  Xg,  x^,...xN  in  D( A)  such  that 


Vi  =  1 , . .  .  ,N  ,  x^  -  U(hi)xi_1  »  h^^ 


where  hi  =  t^  -  t^_.j.  We  then  define  v 


°,z,xr 


(,5)  Vo.z.x0(0)  ■  -O'  Vl  -  ’ . .  Vt  e  1*1-1'  'i1'  ''».x.x0lt)  •  *1  • 

Remark .  If  z  =  0,  the  relation  (14)  always  defines  a  sequence  (x^)  such 
that  =  U(h^ )U(h£_i ) , . . . ,U(hj )xg  which  are  the  expressions  we  are  mainly 
interested  in.  However,  many  situations  involve  perturbed  schemes  like  (14) 
where  z  is  not  zero  but  of  small  norm  in  L^OjT;  X)  and  our  method  applies 


as  well  to  this  more  general  situation. 


Our  goal  is  to  describe  the  behavior  of  v 


when 


n  n  n 
0  ,z  ,x 


lim  | on |  =  lim  |zn(t)|dt  =  lim  |xn  -  xQ 
n+®  n-*®  n+® 


=  0. 


We  expect  it  to  tend  to  the  solution  (if  it  exists)  of 


(  16) 


—  +  Au  =  0  on  (0,T),  u( 0 )  =  x.  e  D( A) 
at  U 


By  solution  of  (16),  we  mean  a  function  u  e  C([0,T);  X)  such  that  there 
exists  a  sequence  of  step-functions  uP  (0,T]  -*■  X  satisfying 


(17) 


Vj  =  1,...,NP,  Vs  e  ]SP_1f  SP],  u1 


(  s )  =  up 


where 


( 18) 


0  =  s0  <  s?  <  ‘ 


.  <  sp  <  T  <  sp  ,  =  s?  -  s?  . 

NP-1  NP  3  3  3-1 


(19) 


P  P 
ur  -  u . 


2 - tl  +  Aup  = 


iP 

3 


wP 

j 


(u^  =  up(0)) 


(20 


N 

)  0  =  lim  sup  |uP(s)  -u(s)|  =  lim  \  Jl?|wP|  =  lim  max  &P 

P+®  0<s<T  p-K”  j=1  3  3  p-w  j^jjP  3 


We  will  refer  to  such  a  solution  as  a  mild  solution  of  (16)  (see  [1]). 

Theorem  1,  Let  u  be  a  mild  solution  of  (16).  Let  v11  =  v  be 

-  n  n  n 

a  ,z  ,x 

sequence  of  functions  defined  by  (12)  -  (15)  where 


lim 

n+® 


lim  / 
n+® 


l  |zn(t)|dt 


lim 

n+® 


x  -  x. 


=  0 


Then,  if  (9),  (10),  (11)  hold 


lim  sup  |u(t)  -  v  ( t ) |  =  0 
n*®  0<t<T 


-4- 


,»■>  .v.  - . 


.  "  %V-V  . V  :  . 


.  ;  *  *  .4  *\  V  ,  •*“-**’.  •*  ~ 

>-1  2  .  _ _  t-  l.  H  w  w  .  «Ll. 


In  the  case  when  the  range  of  I  +  Aa  contains  D(A)  for  all  X  >  0 


thanks  to  Crandall  -  Liggett  [5] ,  one  knows  that 


(21)  Vx  e  D ( A )  ,  Vt  >  0,  S(t)x  =  lim  (i  +  -^A )  nx  exists, 

n*»  a 

and  S(t)  is  a  semigroup  of  nonlinear  contractions  from  D(A)  into  itself. 
Moreover,  by  definition,  t  S(t)Xg  is  a  mild  solution  of  (16)  for  all 
T  >  0  . 


In  that  situation,  as  an  immediate  corollary  of  Theorem  1,  we  have 
Proposition  1.  Assume  A  satisfies  (9),  (21)  and  that  (10),  (11)  hold.  Then 
is  a  sequence  such  that 


if  ( hn  J 

1  1<i<N° 


(22) 


lim  (max  h. )  =  0,  lim  £  h.  =  t 


n>"  l<i<Nn 


n->-oo  i=  1 


we  have  for  all  xQ  e  D( A) 


(23) 


lim  u(hnJ  u(hnn-1)  ...  u(hI|)x0  -  S(t)xc 


and  the  convergence  is  uniform  in  t  on  bounded  intervals. 

The  proof  of  Theorem  1  is  based  on  the  next  lemma  which  is  an  adaptation 
of  technics  used  in  [6]  and  [8]. 

Lemma  1 .  Let  up:  [0,T]  +  X  satisfying  (17)  -  (19).  Then,  for  all  e  >  0 
there  exists  n  =  n  (e,p)  >  0  such  that  for  all  functions  v 

°'z'x0 

satisfying  (12)  -  (15)  with  |a|  <  r\  we  have 


Vy  e  D(A),  Vi  =  Vj  «  0,...,nP 


x . 

X 


vi 


(24  ) 


+  l  \[  |zj  +  e)  +  l  *Pk|w£|  , 


-6- 


where  a  =  max  *  ■  and 

k3<np  3 


(25)  M  ( y )  =  max  { ) Ay | ,  max  h  1 |y  -  U(h  )y|}  . 

1<i<N 

Proof .  Lee  us  denote  Ah  =  h_1(I-U(h)).  By  (11),  for  t  >  0  and  p  fixed 
there  exists  0  =  ri(e,p)  >  0  such  that 


h  e  ] 0 , n]  =>  Vj  =  1,...,NP,  |Ahup  -  Aup|  <  E 


Now,  we  will  assume  |o|  <  r)  and,  since  p  is  fixed  we  will  drop  the 


indexation  by  p. 


Let  us  denote  a.  .  = 


i,j  "  K  "  ujl  and 


let  us  establish  (24)  in  four 


steps.  We  fix  y  e  D(A) 


i)  estimate  of  aQ  j since  A  is  accretive  and  by  (19) 

|Uj  -  y|  <  |u^  -  y  +  ^(Au^.  -  Ay)  |  <  |uj_1  -  y|  +  ^ ^ | Ay |  . 


By  induction,  for  all  j  >  1 


|uj  -  y|  <  |uQ  -  y |  +  ^ikwk  +  s j  | Ay |  . 


Hence,  for  all  j  >  0 


(27)  a0, j  <  lx0  -  yl  +  luo  ■  yl  +  sjlAyl  +J1Akwk 

which  implies  (24)  for  i  =  0  and  j  =  0,...,NP  . 
li)  estimate  of  ax,Q: 

|xx  -  y |  <  |Xl  -  U(hi)x._1|  +  |u(hi)x±_1  -  U(h.  )y|  +  |u(h.)y  -  y| 
<  h±  | z±  |  +  |xi_1  -  y|  +  h.M0(y)  . 


By  induction,  for  all  i  >  1 


Ixi  -  vl  <  lx0  -  y|  +  l  >\lzkl  +  tiMo(y) 

k=  1 


Hence  for  all  i  >  0 


(28)  ai/0  <  lu0  '  yl  +  lx0  “  y|  +  l  hklzkl  +  tiM0(y)  . 

k=  1 

iii)  estimate  of  ai , j  in  terms  of  a^_1fj  and  a^  j-i:  for  « 

i,  j  >  1 

(29)  [l  +  £  h.^ja.  .  =  |x  -  u.  +  i  .h.  Vx.  -  U(h.)u.  +  U(h.)u.  -  u.i 

v  3x1,3  i  3  ji'i  i  j  x  j  j 


Using  (14),  we  have 


lx.  -  U(h . )u .  I  <  a .  ,  .  +  h,  I z .  |  . 

1  l  l  j  1  1-1,3  i  1  i  1 

Plugging  this  inequality  in  (29)  together  with  (26)  applied  with  h  =  h^  <  n 
leads  to 

(30)  [1  +  £.h-1 ja.  .  <  £.h-1  a.  ,  +  i. IzJ  +  £.e  +  lx,  -  u.  -  £.Au. I  . 

v  j  i  i.-  3  ^  3  j  jj 

Now  by  ( 19 ) 


lxi  -  uj  -  Vv  <  ai,j-i  +  Slwjl  • 


Finally,  (30)  and  (31)  yield 

£  h  h 

<33>  Vj  'irrrn-i.j  'drai,j-i  *irrr-  *  |wj I  +e)- 

i  3  i  3  13 

iv)  induction  procedure:  we  assume  that  (24)  holds  with  (i,j)  replaced 
by  ( i— 1 , j )  and  (i,j).  Plugging  it  into  (33)  gives  (24)  after  a  easy  computa¬ 
tion.  Thanks  to  the  two  first  steps,  we  deduce  that  (24)  holds  for  all 
(  i ,3  )  • 


Proof  of  Theorem  1.  Let  v  =  v  be  a  sequence  of  functions  defined 

-  n  n  n 

0  ,z  ,x 

by  (12)  -  (15)  where 


-7- 


lim  | on |  =  lim  |zn(t)|dt  =  lim  |xn  -  xQ|  =  0  . 

n+<»  n+®  n+®  ^ 

Let  u  be  a  mild  solution  of  (16)  and  let  u^1  be  an  approximated  sequence  as 
in  (17)  -  (20).  Let  us  denote 


By  lemma 


|op|  =  max  £p  . 

K:<NP  : 

1,  for  all  p  >  1  and  e  >  0,  there  exists 


Vn  >  N(p)  V  1<i<Nn  V  1<j<NP  Vy  e  D( A ) 


<  xn  -  y|  +  |up  -  y|  +  {(t"  -  sp)2  +  |an|t"  + 


N ( p )  such  that 


where 


Using 


M^ty)  =  sup  | h  \y  -  U(h)y)  j  . 

he(0,1) 

|vn(t)  -  u ( t ) |  <  | v n ( t )  -  uP( t) |  +  |uP(t)  -  u( t ) | 


we  deduce  that  for  all  p  >  0  and  e  >  0: 


lim  sup  (sup  | 

n+<*.  o<t<T 


'  '  I  / 


0<t<T 


-  y| 


1%  -  y  I 


|ap|  +  t|op| 


V2Moo(y)  +  eT  +  £  *p 
k=1  * 


1^1 


Now  letting  p  tend  to  00 ,  then  e  tend  to  0  yield  (according  to  (20)) 


lim  sup  (sup  |vn(t)  -  u ( t ) | )  <  2 | x  -  y)  . 
n-*-®  0<t<T 


Since  xQ  e  0(A)  and  y  is  arbitrary  in  P(A),  the  conclusion  of  the  theorem 


follows . 


3.  The  counterexample. 

Let  X  =  Cq(R)  the  space  of  continuous  functions  u  on  R  such  that 

lim  u(x)  =  0,  equipped  with  the  norm 

x>±® 

II  ull  =  sup  |  u(  x)  |  . 

xeR 


we  define 


D( A)  =  {u  e  C0(R);  u1  e  CQ(R) } 


Vu  e  D(A),  Au  =  u’ 


Then,  A  is  the  infinitesimal  generator  of  the  semigroup  S(t)  on  X 


defined  by 


Vu  e  Cfl ( R) ,  S(t)u(x)  =  u(x  -  t) 


which  means  that  one  has 


Vu  e  D ( A )  ,  1  im  —  *-■ - S(t)u(x)  _  ui(xj  uniformly  on  R. 

t  +  0 

Now,  let  un  e  X  and  a;  (0,1)  +  R  such  that 


lim  a(h)  =  0  . 
h+O 


We  set 


Vh  e  (0,1),  Vu  e  X,  U(h)u  =  S(h)u  +  a(h)(u0  -  s(h)u0) 


Then  we  have 


Theorem  2. 


l)  Under  assumption  (34) 


Vh  e  (0,1),  Vu,u  II U  (  h  )  u  -  U  (  h  )  u  II  <  Hu  -  ull 


VX  >  0,  Vu  e  X  limfl  +  I-U(h) ) j-1u  =  (I  +  XA)-1u 
h  +  o'  h 


-9- 


ii)  One  can  shoose  u0  e  X,  a  satisfying  (34)  and  a  sequence 


{h”,...,hn  }  such  that 
1  n 


)  hn  =  T,  lim  max  hn  =  0 


n-*00  1<i<N 


lim  ilu(hn  )  u(hn  ) 
n+®  Nn  Nn-1 


•  =  +« 


Proof .  The  first  statement  of  i)  is  immediate  from  the  definition  (35)  and 
the  fact  that  S(h)  is  nonexpansive  in  X. 

Let  u  e  X,  and  uh  =  [I  +  Ah-1 (I  -  U(h))]-1u,  that  is 

(36)  u^  +  Ah  1  ( u,  -  U(h)u^)  =  u 


or,  by  (35) 


+  Ah  1(uh  -  SOOu^  -  a(h)(uQ  -  S(h)uQ))  =  u 


which  can  be  written  as 


(37)  u^  -  a(h)uQ  +  Ah  1  (u^  -  a(h)uQ  -  SfhHu^  -  o(h)u0))  =  u  -  a(h)uQ  . 

Setting  =  [I  +  Ah  '(I  -  S(h))]  \  (36)  is  equivalent  to 

(38)  uh  =  o(h)uQ  +  J^(u  -  a(h)uQ)  . 

Since  S(h)  is  a  linear  semigroup  of  contractions,  is  also  nonexpansive 

and  satisfies  (see  for  instance  (9]) 


lim  J^u  =  (I  +  Aa)  ^u  . 
h+0 


But  now,  since  is  nonexpansive,  by  (38) 


Vo  I  <  2a(h)KI  * 


This  together  with  (39)  and  (34)  proves  the  assertion  i). 

Now  we  choose  a  sequence  of  numbers  (h?)  in  (0,1)  satisfying 

1<i<N 

Nn 

(40)  l  hn  =  T  >  0 

i=1  1 


(41 ) 


Vn,  i  /  j  =>  h£  /  hj 


(42) 


—2—  <  min  hn  <  max  h"  <  —  . 
n+1  1<i<Nn  1  1<i<Nn  1  n 


The  conditions  (41)  -  (42)  insure  that  all  the  numbers 

{h^,  1  <  i  <  N11,  n  >  1}  are  different.  This  will  allow  us  to  assign  the 

values  ot(h")  as  we  wish.  We  do  it  as  follows, 
i 

Set  t  “  1  h^,  tQ  =  0,  and  choose  uQ  e  Cg(R)  such  that 

k=0 

Nn 


(43) 


i=1 


tend  to  ®  when  n  tends  to  «.  Actually  any  uQ  e  Cq(R)  whose  total 
variation  on  [0,T]  is  infinite  will  be  convenient.  Indeed  0n  is  the  total 
variation  of  the  step  function  Ug  defined  by 


Vi  =  1,...,Nn,  Vt  e  ]t”_1#  t^J,  u£(t)  =  Ug(t”),  Ug ( 0 )  =  uQ ( 0 ) , 

and  Ug  converges  uniformly  on  [0,T]  to  Ug  when  n  tends  to  °». 

Therefore  the  total  variation  of  Ug  cannot  stay  bounded  if  the  variation 
of  Ug  is  infinite. 

Now,  for  reasons  that  will  become  clear  later  we  choose  a:  (0,1)  ♦  R  as 


follows: 


if  h  =  h£,  a(h)  =  9"  7  sign  (uQ(tJ)  -  «0(tJ_t) ) 
if  h  IS  {h?;  1  <  i  <  Nn,  n  >  1},  a(h)  =  0  . 


Note  that  a  satisfies  (34). 

We  claim  that  with  this  choice,  if  we  denote 

fn  =  u(hn  J  u(hn  )  ...  u(h"lo, 

vn'  v  n  '  v  1 ' 

N  N  -1 

then  fn(T)  tends  to  00  when  n  tends  to  “  and  therefore 

(45)  lim  U  f nU  =  “  . 

00 

n+» 

Indeed,  if  f£  =  u(h£)  u(h”  ^  ...  u(h^)o  for  k  >  1  and  = 

definition  (35)  we  have 


.  .(hy  (u0  -  s(hy  u0) 


By  induction,  we  check  that 


f"  =  l  ot(hn)  [S(t"  -  tn)u 

k  i  k  i  0 

i=1 


s(tk-  Vi’V 


so  that 


fn  =  i  o(h")[S(T  -  tj)u  -  S(T  -  tn  • 

i=1 

But  S ( T  -  t")uQ(x)  =  Uq(x  -  T  +  t?),  so  that  according  to  (44) 

(46)  fn(T>  =  9^1/2  l  |u0(t")  -  u0(t^_1)|  =  e^2. 

k=  1 

From  (43),  (46)  we  obtain: 


lim  fn(T)  =  lim  0  ® 


Remarks . 


i)  The  same  negative  conclusion  could  be  reached  by  choosing 

U(h)u  =  (I  +  hA)  ^  +  a(.,)(u  -  (I  +  hA)  ^u^) 

which  is  nothing  but  the  resolvents  of  the  operator  u  A(u  -  a(h)Ug).  Note 
that  U(h)  as  defined  by  (35)  is  the  value  at  h  of  the  semigroup  generated 
by  the  operator  u  +  A(u  -  a(h)Ug). 

ii)  According  to  this  counterexample,  the  conclusion  (23)  dramatically 
fails  as  soon  one  switches  from  a  regular  step-size  to  a  very  slightly 
irregular  path.  In  view  of  (41),  (42)  we  need  all  the  h?  to  be  different 
but  they  can  be  as  close  to  each  other  as  we  pleased  (see  (41),  (42)). 

iii)  It  would  be  interesting  to  know  whether  there  exists  a  condition 
stronger  than  (4)  but  weaker  than  (2)  such  that  the  expected  conclusion  (8) 
holds.  Let  us  mention  some  positive  results  with  variable  step-size  obtained 
in  [10]  under  stability  assumptions  on  the  mappings  U(h). 

iv)  If  B,  C  and  B  +  C  are  infinitesimal  generators  of  linear 
semigroups  of  contractions  and  if  one  sets 


A  =  B  +  C 

U(h)  =  ( I+h  B)-1  ( I+h  C)_1  or  U(h)  =  SB(h)  SC(h), 

B+C 

one  easily  verifies  that  assumption  (11)  holds.  Since  S  (t)  uQ  is  a  mild 
solution  of  ^7  +  Au  =  0,  u(0)  =  u^  for  all  Ug  (see  [1]  ),  for  both  choice 

one  has,  as  a  consequence  of  Theorem  2.1 


lim  u(hn  )  u(hn  )  ...  U^h1?)  u.  =  SB  C(t)  un  . 

'  n  n  .  v  1 '  0  0 

n+“  N  N  -1 

v)  The  same  conclusion  holds  if  A  =  B  +  C  where  B  is  the 
infinitesimal  gnerator  of  a  linear  semigroup  of  contractions  and  c  is  a 


continuous  accretive  operator  defined  on  the  whole  space. 


References 


[ 1 }  Benilan  Ph . ,  M.  G.  Crandall,  A.  Pazy,  Evolution  equations  governed  by 
m-accretive  operators,  to  appear . 

12)  Brezis  H.,  A.  Pazy,  Convergence  and  approximation  of  semigroups  of 

nonlinear  operators  in  Banach  Spaces,  J.  Funct.  Anal.  9  (1972),  63-74. 

[3]  Chernoff  P.  R. ,  Note  on  product  formulas  for  operator  semigroups,  J. 
Funct.  Anal.  2  (1968),  238-242. 

[41  Chernoff  P.  R. ,  Product  formulas,  nonlinear  semi-groups,  and  addition  of 
unbounded  operators,  Mem.  Amer.  Math.  Soc.  140  (1974). 

[5J  Crandall  M.  G. ,  T.  M.  Liggett,  Generation  of  semigroups  of  nonlinear 
transformation  on  general  Banach  spaces,  Amer.  J.  Math.  93  (1971). 

[6]  Kobayashi  Y.,  Difference  approximation  of  Cauchy  problems  for  quasi- 
dissipative  operators  and  generation  of  nonlinear  semi-groups,  J.  Math. 
Soc.  Japan,  27  (1975),  640-665. 

[7]  Kurtz  T.  G. ,  Extensions  of  Trotter's  operator  semigroup  approximation 
theorem,  J.  of  Functional  Ana.,  3  (1969),  354-375. 

[8]  Miyadera  I.,  Kobayashi  Y.,  Convergence  and  approximation  of  nonlinear 
semi-groups,  in  Functional  Analysis  and  Numerical  Analysis,  Japan-France 
Seminar,  Tokyo  and  Kyoto,  (1976);  H.  Fujila  (Ed.):  Japan  Society  for  the 
Promotion  of  Science  (1978). 

[9]  Pazy  A.,  Semigroups  of  linear  operators  and  applications  to  partial 
differential  equations.  Applied  Mathematical  Sciences  44,  Springer- 
Verlag,  New  York  (1983). 

[10]  Schechter  E. ,  Stability  conditions  for  nonlinear  products  and  semi¬ 
groups,  Pacific  S.  of  Math.  85  (1979),  179-200. 

[11]  Trotter,  H.  F.,  Approximation  of  semigroups  of  operators,  Pac.  J.  Math., 
8  (1958),  887-919. 


SECURITY  CLASSIFICATION  OF  THIS  PACE  (TWian  Oil*  Bntorot o 


REPORT  DOCUMENTATION  PAGE 


I.  REPORT  NUMBER 


4  TITLE  (end  Submit)  S.  TYPE  of  REPORT  a  period  covered 

Summary  Report  -  no  specific 

About  Product  Formulas  With  Variable  Step-Size  reporting  period 

s.  performing  org.  report  number 


READ  INSTRUCTIONS 
BEPORE  COMPLETING  PORM 


QKMENT'S  CATALOG  NUMBER 


I.  REQPWEN  ~  S  I 


7.  AUTHOR/ a.) 


Michel  Pierre  and  Mounir  Rihani 


»•  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 


Mathematics  Research  Center,  University  of 
610  Walnut  Street  Wisconsin 

Madison,  Wisconsin  53706 


I  I.  CONTROLLING  OFFICE  NAME  AND  AOORESS 

U.  S.  Army  Research  Office 
P.O.  Box  12211 

Research  Triangle  Park,  North  Carolina  27709 


.  MONITORING  AGENCY  NAME  A  ADDRESS/ff  df  Karan  f  from  Controlling  Ollleo) 


I.  CONTRACT  OR  GRANT  NUMBER^ 


DAAG29-80-C-00  41 


<0.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 

Work  Unit  Number  1  - 
Applied  Analysis 


12.  report  oatb 

January  1985 


IS.  NUMBER  OF  PAGES 

14 


IS.  SECURITY  CLASS,  (ot  Ihto  rmpo rt) 

UNCLASSIFIED 


ISa.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


IS.  DISTRIBUTION  STATEMENT  (ot  (Ala  Roport) 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ot  tho  obotrael  an  farad  In  Block  30,  It  dtttoro nf  from  Roport) 


19.  KEY  WORDS  ( Contlnu e  on  r»r«r««  aide  it  necemamry  and  Identify  by  block  number) 

Product  formulas,  nonlinear  semigroups,  accretive  operators, 
approximation  scheme,  stability  and  convergence  » 


20.  ABSTRACT  (Continue  on  reveree  aide  tt  naceaamry  and  identity  by  block  number) 

A  continuous  semigroup  of  linear  contractions  (S(t))  t  >  0  on  a  Banach 
space  X  can  be  generated  through  the  product  formula 

(1)  S(t)x  -  lim  uf  — ) nx  Vx  e  X  Vt  >  0 

n 


DD  |  j»H^j  1473  EDITION  OF  t  NOVAS  IS  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  /Whan  Dole  Bntorod) 


20.  Abstract  (cont.) 


where  (u(h))h>0  is  any  family  of  contractions  on  X  whose  right-derivative 
at  h  =  0  coincides  with  the  infinitesimal  generator  A  of  S(t),  that  is 
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h+0  h  t*0  11 


Actually  formula  (1)  extends  to  the  case  when  S(t)  is  a  semigroup  of 
nonlinear  contractions  "generated"  by  an  m-accretive  operator  A  and 
condition  (2)  may  also  be  weakened  to 


(3) 


i  -  u(h)  hiy  A 
h 


in  the  sense  of  graphs. 


We  look  here  at  the  same  formula  (1)  when  the  regular  step-size  t/n  is 
replaced  by  a  variable  step-size,  namely 

(4)  lim  u(hn  )  u(hn  ]  ...  u(h?)x  =  S(t)x  Vx  e  D( A) 

n+o°  Nn  Nn-1 

where  )  hn  =  t  and  lim  f  max  hnl  =  0. 
u  n  1  v  n  i  * 

14i<N  nx»  Ki<N 

Surprisingly,  it  turns  out  that  (4)  fails  under  assumption  (3);  we 
exhibit  a  counterexample.  However,  we  prove  that  (4)  holds  true  under  the 
stonger  assumption  (2). 
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